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We evaluate the dispersion relation for massless fermions, described by the Dirac equation, and 
for zero-spin bosons, described by the Klein-Gordon equation, moving in two dimensions and in the 
presence of a one-dimensional periodic potential. For massless fermions the dispersion relation shows 
a zero gap for carriers with zero momentum in the direction parallel to the barriers in agreement 
with the well-known "Klein paradox". Numerical results for the energy spectrum and the density 
of states are presented. Those for fermions are appropriate to graphene in which carriers behave 
relativistically with the "light speed" replaced by the Fermi velocity. In addition, we evaluate the 
transmission through a finite number of barriers for fermions and zero-spin bosons and relate it with 
that through a superlattice. 
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I. INTRODUCTION 

The recent realization of stable single layer and bilayer 
carbon crystals (graphene) has aroused considerable in- 
terest in the study of their electronic properties^ ^ . These 
materials have unusual properties that may allow the de- 
velopment of carbon-based nanoelectronic devices. The 
behavior of charge carriers in wide single-layer graphene 
sheets is that of chiral, "relativistic" massless particles 
with a "light speed" equal to the Fermi velocity of the 
crystal {vp ^ c/300) and a gapless linear dispersion close 
to the K and K' points^ ^ ^ . One consequence of that is 
that single- layer graphene displays an unusual quantum 
Hall effect, in which the quantum Hall plateaus occuil^at 
half-integer multiples oiAe^/h. The absence of a gap and 
the chiral nature of the electronic states, in both single- 
layer and bilayer graphene, is at the root of phenomena 
such as the Klein tunneling^ ^ which is the perfect trans- 
mission of carriers, upon normal incidence, through a 
potential barrier. The study of this effect is relevant to 
the development of future graphene-based device^. 

From the standpoint of basic research, due to its lower 
" light speed" graphene can also fulfill the role of a testbed 
for the investigation of relativistic quantum effects, such 
as Zitterbewegung^^ and pair creation^^. Moreover, the 
role of chirality in the electronic properties of massless 
fermions can be assessed or further appreciated by con- 
trasting the behavior of chiral, massless fermions with 
that of non-chiral, massless zero-spin bosons. 

In this work we study the dispersion relation for two- 
dimensional fermions, described by the Dirac equation in 
the presence of a one-dimensional (ID) periodic poten- 
tial. This can be realized in a periodically gated graphene 



layer. In doing so we extend the treatment of Ref. |T3] 
which considered motion only along one direction. For 
a clearer understanding we contrast these results with 
those for bosons of zero spin that are described by the 
Klein- Gordon equatiorP^l. Here too we include the mo- 
tion parallel to the barriers, which, to our knowledge, has 
not yet been reported. In addition, we briefly present re- 
sults, in Sec. II, for the transmission through a single 
barrier. In Sec. HI we present the dispersion relation 
and transmission for either superlattice case and relate 
the transmission to that through a single barrier. A sum- 
mary and concluding remarks follow in Sec. VI. 



II. TUNNELING THROUGH A SINGLE 
BARRIER 

A. Bosons 
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FIG. 1: (Color online) ID Potential barrier V(x) of height V 
and width W 



We first calculate the transmission of zero-spin bosons 
through a single barrier of height V and width W, shown 
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in Fig.[T] using the Klein-Gordon equation in two dimen- 
sions 



1 



[{E -V{x)f -m^c^]ilj{x,y). (1) 



Since the hamiltonian H commutes with py we look for 
solutions in the form ip{x^y) = ip{x)e^^yy . With this 
substitution the resulting equation for iIj{x) is solved by 
il){x) = e^*^^^ outside the barrier with kx = E' /he and 
by il){x) = e^*^^ inside it with 



n^f^ ^^^^^^ ^^^^ ^ ^ ^ ' (2) 

and with E^ given by 

_ ^2 _ ^2^2^2 _ ^2^4 ^ ^2^2^2^ (3) 

Then following the standard procedure of matching the 
solutions and their derivatives at the interfaces of the 
regions I and II, and II and III, shown in Fig. [l] leads to 
the transmission 
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(4) 



If K becomes imaginary, then K is replaced by i\K\ and 
sin(i^W) by zsinh(i^iy) in this equation. Notice that 
in contrast to the non-relativistic case, K depends on 
ky. The result for the non-relativistic case is obtained by 
inserting ky = 0, k^ = ^ImEjY?^!^, and K = \lm{E - 
V)/h^Y^^ in Eq. (4). 

In Fig. [2] the transmission through a single potential 
barrier is plotted for zero-spin bosons with the same ve- 
locity (c = lO^m/s) and parameters for the potential 
barrier as for the massless fermions in graphene. In do- 
ing so we will be able to better compare the transmission 
for bosons with that for fermions. As shown in Fig. [4j in 
which we plot slices of this plot at constant ky^ the trans- 
mission is perfect for certain values of kx although the bo- 
son has an energy that is much lower than the height of 
the potential barrier V (if ky = then E = V = 50meV 
for kx = 0.076nm~^). Notice that these transmission res- 
onances below the barrier height do not occur in the non- 
relativistic case. Moreover, the transmission depends on 
the wavevector component parallel to the barrier, i.e., on 
ky since K does, cf. Eq. (2), whereas it does not in the 
non-relativistic case. 

Since K is always real for m = 0, the sine term in 
Eq. (4) only introduces a modulation of the transmission. 
The minima in the transmission arise from the (/c^ — 
K'^)/2kxK term. This is because at the vicinity of £^ = 
the solutions inside the barrier approach those of a 
free particle with zero energy and thus with /c = 0, i.e. 
infinite wavelength. Together with the continuity of ^ 
and of its derivative, this acts to reduce the amplitude 
of the transmitted wave. For the non-relativistic case 
and E < V ^ K becomes imaginary and the sine-term 
becomes a hyperbolic sine, which can become quite large 
and reduce significantly the total transmission. 
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FIG. 2: (Color online) Contour plot of the transmission T of 
relativistic bosons through a barrier as a function of kx and 
ky for y = 50 meV, c 10^ m/s, VI/" = 50 nm, and m = 0. 



B. Fermions 

In two dimensions spin- 1/2 fermions are described by 
the 2D Dirac equation: 



[c((j • p) + mc^cFz] = {E - V)^, 



(5) 



where a = {(ix^o-y) and are the Pauli-matrices. In- 
stead of imposing a certain form of solution, without giv- 
ing any details, to Eq. (S)*", we look for solutions in the 
form ip{x^y) = ip{x)e^^yy ^ where the wave function ip{x) 
is a two-component spinor ?/^(x) = [^/^^^(x), ?/^/(x)]-^. Sub- 
stituting this form of solution and setting E' = E — V{x) 
gives 



dilJi{x)/dx + kyijji{x) = —{E' - mc^)\l)u{x), 
di\)u{x)ldx - kyi\)u{x) = Yc{E' + mc^)ilJi{x). 



(6) 



These coupled equations admit solutions of the form 

ij{x) = (""A e^^-^ (7) 



where kx = {l/hc)[E'^ — h^kyC^ — w?c^Y^'^ must be re- 
placed by = (l/hc)[{E-VY-h^klc^-m^c^Y^'^ inside 
the barrier. Then substituting this form in Eq. ([5| gives 



E' — mc^ —hckx- 
—hckx,^ E' + mc^ 



0, 



(8) 



where kx^± = kx :^iky. A non-trivial solution demands a 
zero determinant and this leads to 



2^4 



E'^ = h^c^ki + h^c^k^y + m^c 



(9) 
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and to 



hck. 



Ul, 



E' >0 



E' < 0. 



(10) 



As detailed in the Appendix, the general solution is a 
linear combination of two independent solutions (j)a{x) 
and (l)b{x) of this form. In contrast with the boson case 
of Sec. II. A, the continuity of the derivative of the spinor 
components for Dirac fermions is not necessary. Thus, 
matching only the solutions at the interfaces between the 
regions I and II and II and III gives, for £^ > the 
transmission as 



T{kx', ky) 



n -1 



1 sin2(KW) 



(11) 
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FIG. 3: (Color online) Contour plot of transmission T of 
massless relativist ic electrons through a barrier with m = 0, 
y = 50 meV, c = 10^ m/s, W = 50nm 
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with d = {\E\ + mc^)/hc and D = (\E -V\^ mc^)/hc. 
Similar expressions are obtained for < E < V and 
E <0, d. Appendix. 

For \E — < h^kyC^ + m?c^, K is imaginary and 
the transmission is given by Eq. ( pTj ) with K replaced 
by iK' = i\K\ and sin(i^rW) by zsinh(i^W). A contour 
plot of the transmission is given in Fig.[3j As can be seen, 
for normal incidence, i.e., for ky = 0, the transmission is 
equal to 1. Contrasting Fig. [3] with Fig. [2] for bosons 
we see a significant gap in Fig. |2] for normal incidence 
that is not present in Fig. [3j Another way to contrast 
the results of the two cases is shown in Fig. [4] where 
slices of the transmission across the planes ky = and 
ky = 0.05 nm~^ are shown for bosons (green curves) 
and electrons (red and magenta curves). The perfect 
transmission (T(kx^ ky) = 1) for normal incidence results 
from the chirality of the carriers which prevents them 
from being refiected by a barrier, see also Ref.|7l Because 
massless bosons are not chiral, they are not subject to 
such a constraint. 
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FIG. 4: (Color online) Plot of slices of the transmission co- 
efficient for bosons (green curves) and electrons (solid red 
and dashed magenta curves) taken at constant ky = and 
ky = 0.05 nm~^ from Figs. [5] and [s] respectively. 



III. SUPERLATTICES 

We now consider ID superlattices (SL), that are peri- 
odic structures consisting of N identical units of length 
L. Each unit consists of a rectangular barrier and a rect- 
angular well whose widths are denoted, respectively, by a 
and b. We first consider fermions described by the Dirac 
equation, Eq. ^ . For motion of electrons along only one 
direction (x), such a calculation was done by McKellar et 
al^^. We generalize this earlier work to two dimensions. 

We start with Eqs. taking c ^ lO^m/ for 

our numerical results, a value that is appropriate for 



graphene. Equation ^ can be written as 

^(-) = (til) • 

Using Eq. ^ we have two independent solutions (j)a{x) 
and (pbix) in each unit cell. The solution in the nth unit 
cell n, written as tpn{x) = V^(x + (n — 1)/), is a linear 
combination of these two 

*'-<^> = (til) + (til) ■ !»' '1- ("> 
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It is convenient to define the matrix 



(15) 



Then the continuity of the wave function leads to the fol- 
lowing relation between the coefficients of the unit cells: 



(16) 

(17) 



and the transfer matrix T is given by T = ll(0)~^ri(/). 

It can be proven^^ that det[T] = 1. Assuming periodic 
boundary conditions we can write 



V^(x + L) = V^(x), 

f CtN-\-n\ I Ct 



and obtain 



\bN-\-n J 



1. 



(18) 
(19) 

(20) 
(21) 



It follows that the eigenvalues of T are e^^*^/^, with n an 
integer, and its determinant is equal to 1. The dispersion 
relation is then given by the trace of T as 



2cos{kJ) =Tr(T), 



(22) 



where kx = 27rn/L. Notice that this kx^ or Bloch 
wavevector, expresses the periodicity of the structure and 
should not be confused with the kx of Sec. 11. 



A. ID Superlattices in graphene 

The dispersion relation is given by Eq. (22) and the 
transfer matrix T by 



(23) 



where Qk and fl^ are the matrices in the barriers and 
wells, respectively. With the help of the Appendix we 
can specify the matrices fix and ft^. For E > V the 
dispersion relation becomes {k, = E' /he) 

2 cos{kxl) = 2 cos(i^a) cos{f<ib) — G sin(i^a) sin(/^6) , (24) 

where 

kK kK' ^ ' 



G = 



dp 



By setting fcj, = in Eqs. (24) and (25) we recover the 
result of Ref. [13 for ID motion of Dirac particles. For 



{E - Vf < m^d^ + fi.2fc2g2 and E - y > we have an 
imaginary K and simply make the substitution K — )• i\K\ 
in Eq. (24). For E < the result for G is the same 
whereas for F > £J > it becomes 



G 



dP 



■fc2)(«2. 



KndD 



2fc| 
kK' 



(26) 



The first two minibands of the dispersion-relation for a 
ID graphene-based SL are plotted in Fig. [5] Notice that 
the form of the dispersion relation in the y direction, in 
contrast with the non-relativistic case, depends on the x 
component kx of the wave vector and is not the energy 
of a free particle. As a result, we cannot split the energy 
of the minibands in the manner E^ = En{x) + chky. If 
ky is sufficiently large, this dependence goes away and 
we have again a linear relation between the energy and 
ky, independent of the kx component. Notice also that 
for ky = (and m = 0) the linear dispersion relation 
E{kx,0) = vpfikx is unaltered by the periodic potential 
F, it is only shifted by V/2. This is a consequence of the 
Klein paradox. 

In Fig.|6]we show three slices of the dispersion relation 
of Fig. pj taken at constant ky = 0, 0.066, 0.132 nm~^ 
(red solid, green dashed, and blue dash-dotted curves, 
respectively). We see that the ky component is also gen- 
erating a band gap. 
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FIG. 5: (Color online) Dispersion relation in a ID graphene 
SL. Only the first two bands are shown for a = 6 = 10 nm, 
V = 100 meV, and m = 0. 

It is instructive to take the Kronig- Penney limit of the 
dispersion relations given above. For very thin and high 
barriers, such that a ^ and V ^ oo but with their 
product (5* = aV) kept a constant we have 



K 



E^ — h?c^ky'^ — im?d^ 

he 

^{E - V)2 - h^c^ky^ - 

he 



v_ 
he 



(27) 
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occurs by "inserting" small gaps between regions of per- 
fect transmission. The resonances also show a very slight 
increase in the transmission. This behavior remains es- 
sentially unaltered if the number of units is increased 
and is similar to that reported previously in the context 
of resonant tunneling through a double barrieil^. 
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FIG. 7: (Color online) Contour plot of the transmission T 
of relativistic fermions in a structure consisting of ten units 
with c — 10^ m/s, V = 50 meV, a = 50 nm, b — 100 nm, and 
m = 0. 



FIG. 6: (Color online) Slices of the dispersion relation shown 
in Fig. [5 (a) is for constant ky = 0, 0.066, 0.132 nm~^ for the 
solid red, dashed green, and dash-dotted blue curves, respec- 
tively, and (b) for constant /c^ = 0, tt// for the solid magenta 
and dashed green curves, respectively. 



K: E'^ — hP'c^ky'^ — im?c^ K 
d ~ 1^1 +mc2 " ~ 



Khc 



D 



\E' 



^ 1. 

(28) 

For V ^ 00 only the cases with positron waves in the 
barrier are allowed. Then the dispersion relation becomes 

1) \E\ > mc^ 

S S 
2cos(kxl) = 2 cos(-— ) cos(ai:/) + G+ sin( — ) sin(Ai:/), (29) 

he he 

2) \E\ < me^ 

S S 
2 cos{kxl) = 2 cos( — ) cosh(/^/) + G- sin( — ) sinh(Ai:/), 

(30) 

where G± = hz/d ± d/hz -\- ky/nd. 

To make contact with Sec. II, we now consider the 
transmission through a multi-barrier structure that is a 
good approximation to a SL. In Fig.[7|we show a contour 
plot of the transmission for a structure consisting of ten 
units. Upon contrasting Figs. [3] and [7| it is seen that the 
main effect of having many units is to slightly reduce the 
continuous range of kx for which a perfect transmission 



B. Bosons 

The dispersion relation for a SL and relativistic bosons 
is obtained in a similar way and in full analogy with the 
non-relativistic Schrodinger equation. We use the Klein- 
Gordon equation, Eq. ([T]), and obtain the same dispersion 
relation as in the non-relativistic case but with different 
values for n and namely, 

n = {l/he)[E^ - h^ky - m2c^]^/^ (31) 
K = {l/he)[{E - - h^ky - m^c^]!/^ (32) 



The dispersion relation is 



2cos{kxl) = 2cos{Ka) cos{hzb)- 



^m.{Ka) sin(/^6). 

(33) 

If K becomes imaginary, which happens for {E — V^ < 
m?e^ + h^l^e^ ^ we simply replace K hy i\K\ in Eq. (34). 



In Fig. |8] we plot the first two minibands as a func- 
tion of kx and ky for the parameters shown in the cap- 
tion. The main difference with the non-relativistic case 
is the appearance of a large gap along kx and part of 
ky^ centered around /c^^ = 0. Notice the turning point of 
the band, which leads to an infinite group velocity Vky- 
Such an unusual result is not uncommon in connection 
with the Klein-Gordon equation^^. If we take the mass 
sufficiently different from zero the turning point in the 
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superlattice band structure disappears. Notice the sim- 
ilarity and differences with the relativistic fermion case 
depicted in Fig. [5] The situation is similar between Figs. 
[6] and [6] where the dispersion relation is shown vs kx in 
the left panels and vs ky in the right ones. 
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FIG. 8: (Color online) Dispersion relation for relativistic 
massless zero-spin bosons. Only the first two minibands are 
shown for c — 10^ m/s, V — 100 meV, a = 6 = 10 nm, and 
m = 0. 



Again we make contact with Sec. II by considering 
the transmission through a ten-unit structure described 
above for fermions. A contour plot of the transmission 
Contrasting Figs. [2] and [lO 



is shown in Fig. [T0| Contrastmg l^ igs. and [UJ] we see 
that the main effect of having many units is to reduce the 
continuous range of kx for which a perfect transmission 
occurs, by "inserting" small gaps between regions of per- 
fect transmission, and introduce resonances in the gap 
of Fig. [2] which occurs for 0.5 < kx < 0.014. This was 
not the case for fermions for which the transmission was 
perfect in this range. Another way to contrast the results 
of the two superlattice cases considered above is shown 
in Fig. [TT] where the transmission is plotted vs kx for ten 
units and constant ky. Panel (a) is for normal incidence, 
that is for ky = 0^ the red (straight) curve being for elec- 
trons and the blue one for bosons. Panels (b) and (c) are 
for electrons and bosons, respectively, and in either case 
for ky = 0.02 nm~^. We see again the perfect transmis- 
sion for electrons in panel (a) and the quite different one 
for bosons while for an oblique incidence the difference is 
less pronounced. Notice the quite " square- wave" charac- 
ter of the transmission, which becomes more pronounced 
upon increasing the number of units. Notice also the 
similarity between Figs. |4] and jT Again the differences 
between electrons and bosons are due to the chirality of 
the former. 
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FIG. 9: (Color online) Slices of the dispersion relation for 
relativistic bosons extracted from Fig. [s] (a) is for constant 
ky — 0,0.16 nm~', for the solid red and dashed blue curves, 
respectively, and (b) for constant kx — 0, tt/^ for the solid red 
and dashed green curves, respectively. 



C. Density of states 

An additional way to contrast electrons with bosons 
is to evaluate the density of states (DOS) D{E). In the 
reduced-zone scheme it is given by 



tt/I 



dkx / dkyS{E - En{kx,ky)), 
Jo 

(34) 

where A is the surface area. The integral is evaluated 
numerically by converting it to a sum in the manner 



pTv/l poo ^ C 



where the kx and ky indices take the values 

l---Nx{Ny) 



X y 



(35) 



(36) 



The cutoff C in the ky direction is chosen sufficiently 
large so that there is negligible change in the DOS by 
not considering contributions from ky > C; for the plots 
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FIG. 10: (Color online) Contour plot of the transmission of 
relativistic bosons in a structure consisting of ten units with 
c = 10^ m/s, y = 50 meV, a = 50 nm, 6 = 100 nm, and 
m = 0. 
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FIG. 11: (Color online) Electrons vs bosons. Plots of the 
transmission through a structure consisting of ten units with 
c = 10^ m/s, y = 50 meV, a = 50 nm, 6 = 100 nm, and 
m = 0. 



shown we took C = IOtt//. In addition, we replace the S 
function in Eq. (33) by a gaussian, of width a = 3meV^ 



1 



-[E-Er,iK,ky)f/2a 



\ (37) 



and chose a small but sufficiently large to compensate 
for the discretization of kx and ky. For the DOS of zero- 
spin bosons we divide the result by 4 since we don't have 
the fourfold valley-and-spin degeneracy. We also evalu- 
ate D{E) for a non-relativistic superlattice in the tight- 
binding model 

En{k:r. ky) = + fi'/e^/2m - 2tn cos{kJ), (38) 



where is the middle of the nth miniband and the 
hopping parameter. Notice that we include the free mo- 
tion along the y direction. Then D{E) takes the form 



D{E) = 



^/E -E„ + 2t„ cos{kJ) 



(39) 

and evaluate the integral numerically with tn equal to 
one quarter of the width of the n-th miniband (the pa- 
rameters En and tn we obtained numerically for the same 
shape of superlattice as for the relativistic particles). 

The results are shown in Fig.[l2j The solid red curve is 
for electrons {Dq = AA/tt'^) and the dashed blue one for 
bosons; the left and right scales are different due to the 
absence of spin degeneracy in the latter case. The dotted 
green curve is for non-relativistic electrons in the super- 
lattice in the tight-binding model with Dq = y/2mA/7Th 
and m = 0.067mo. It appears different from the usual ID 
one because we included the motion along the y direction 
but not the one along the z direction^^. For comparison 
we also show the DOS for the Dirac 2DEG in the absence 
of a superlattice by the red dot-dashed curve. The struc- 
ture seen in the DOS is directly related to that of the dis- 
persion relation shown in Figs. |5) [6| [Sj [9] The more pro- 
nounced structure in the case of bosons reflects the wider 
gaps, on the average, at the zone boundaries. In general 
the DOS increases upon approaching the band edges and 
then decreases. Accordingly, there are no peaks in the 
DOS for a Dirac 2DEG, no states with £^ < for a non- 
relativistic particle, etc. 
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FIG. 12: (Color online) Electrons vs bosons. Density of states 
D(E) vs energy E for a superlatice with c = 10^ m/s, 1/ = 100 
meV, a = 5 = 10 nm, m = 0. The red, solid curve is for elec- 
trons and the blue, dashed one for bosons, the latter is shifted 
for clarity. The green dotted curve is for a non-relativistic par- 
ticle (m = 0.067mo) in the tight-binding model, see Eq. (38), 
and the red dot-dashed curve for the massless Dirac electron 
in the absence of a superlattice. 
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IV. SUMMARY 

We obtained the dispersion relation for massless 
fermions and zero-spin bosons using, respectively, the 
Dirac and Klein-Gordon equation with a ID periodic po- 
tential but allowing for motion of these particles in two 
dimensions. In doing so we generalized the treatment of 
Ref. [13] which considered motion only along one direc- 
tion. We also evaluated the transmission through such a 
superlattice and related it to that through a single bar- 
rier. Further, we evaluated the density of states for ei- 
ther superlattice, for a non-relativistic particle, and for 
the Dirac massless fermion. 

The parameters we chose for fermions are appropriate 
to graphene that is currently studied intensively. It is well 
established that electrons in graphene behave effectively 
as chiral, "relativistic" massless (m = 0) particles with 
a "light speed" equal to their Fermi velocity. Though 
the Dirac and Klein-Gordon equations differ significantly, 
we chose the same parameters for zero-spin bosons for 
a more meaningful contrast between the results derived 
from them. The expression for the wavevector Eq. 
(2), indicates that taking m 7^ is equivalent to assum- 
ing a higher effective ky^ through ky^= w?(? jV? + or 
a higher effective mass through m*^ = h^k'^/c^ ^ m?. 
Incidentally, as we pointed out after Eq. (33), taking 
m 7^ gives a finite group velocity Vky in the Klein- 
Gordon case, see Figs. |8] and [9) We notice in passing, 
that despite infinite velocities encountered in problems 
studied with the Klein-Gordon equation, the equation is 
quite versatile as it can also describe waves in plasmas, 
torsion-coupled pendula, pulse propagation in periodic 
dielectric structures as well as tunneling in waveguides of 
photonic band-gap materials, see Ref. [15] and references 
cited therein. 

As expected, we saw that Klein tunnelin^^J^^, the 
perfect transmission of carriers, upon normal incidence 
{ky 7^ 0), occurs not only through a single potential bar- 
rier but also through a superlattice, see Figs. [3] and [TT] 
As pointed out this difference from the non-relativistic 
case occurs because the wavevector K does depend on ky 
in the former case but not in the second. For zero-spin 
bosons this perfect transmission occurs for a highly re- 
stricted range of wavevectors, cf. Figs.[2j[l0j and 11 The 
difference though in the transmission between fermions 
and zero-spin bosons is progressively less pronounced if 
the incidence is oblique {ky 7^ 0), see Figs. [2]|4| 



and 
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Apart from these differences though, it is essential to 
emphasize and realize that a linear spectrum is not an 
imperative requirement in order to have perfect trans- 
mission through potential barriers, for normal incidence, 
i.e., T = 1, whereas chirality is indeed essential, cf. Figs. 
EHanddU 

Finally, we saw the various DOS contrasted in Fig. 12 
Notice that in the DOS for a non-relativistic superlat- 
tice we retained the motion along the y direction, cf. Eq. 
(38). Apart from being meaningful or instructive, the dif- 
ferences between the various cases shown will be reflected 



in the transport coefficients, postponed to a future study, 
and hopefully of use to pertinent experiments. 
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V. APPENDIX 

Below we evaluate the transmission for electrons. The 
solution to Eq. (5) is a linear combination of two inde- 
pendent solutions (t)a{x) and (phi^)- Outside the barrier 
we have 



A+ -A_ '^-^ 



(40) 



where N^^ is a normalization factor, \± and A± are de- 
fined as 



A. 



k±hc 



A4 



K^hc 



(41) 



and Clk^{x) is given by Eq. (15). Inside the barrier 







A+ 


-A_ 


-A. 


_ A+ 







Nkc 



iKoi 



Ke 



E>V 



E<V. 



(42) 



We now use the transfer matrix method and express the 
solutions in terms of (pa and ^5. Referring to Fig. [l] we 
have 



Xx) = Q.kA^){J^Q] 



in region I, 
in region II, 
in region III. 



(43) 



The wave function has to be continuous at x = and 
X = W. This gives 

f)=^k'imK{0)(^'^y (44) 



and 



where T is defined as 



(45) 
(46) 

(47) 



9 



Explicitly we obtain the following results for the matrix 
elements T^n of the 2x2 transfer matrix T. For E > V 
or E < the results are 

Tn = ae""'"^ [|A* + Aj^e"*^^ - |A - Aj^e*^'^] , 

Ti2 = 2iae-''==="'(A - A)(A + A*) sm{KW), (48) 

-'21 = -J-12J -'22 = -til, 

where a = dD/Ak^K. For < i? < V the results are 

Tn = ae""'"^ [-|AA - l^e"*^'^ + |AA* + Ij^e*^*^] , 
Ti2 = 2iae''===^(-A*A* + 1)(A*A + 1) sm{KW), 

rji rT"!* rji rj-f^ 

-1-21 = --I-121 -1-22 = -l-li' 



With reference to Fig. [T] the current is given by 

= 2ecNlXi\Af - \Bf), in region I, 
= 2ecNlX\F\^, in region III. 



(49) 
(50) 



Then the transmission T = T{kx, ky) is 



T 



(51) 



Explicitly, for £^ > V the transmission is given by Eq. 
(11) whereas for < £^ < F it takes the form 



T = [4a^F'^ sm\KW) + cos\KW)] ' 



(52) 



with F' = {K^ + kl){kl + kl)/(PD^ + 1 + 2kl/dD , d = 
{E^mc^)/hc, and D = {\E -V\ ^mc^)/hc. In a similar 
way we obtain the transmission for £^ < 0, which is the 
same as that for > 0, and the other cases listed in Sees. 
II and III. 
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